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We study the physics of four-photon states generated in spontaneous parametric down-conversion 
with a pulsed pump field. In the limit where the coherence time of the photons t^^ is much shorter 
than the duration of the pump pulse At, the four photons can be described as two independent pairs. 
In the opposite limit, the four photons are in a four-particle entangled state. Any intermediate case 
can be characterized by a single parameter x, which is a function of t^^/At. We present a direct 
measurement of x through a simple experimental setup. The full theoretical analysis is also provided. 



I. INTRODUCTION 

Spontaneous parametric down-conversion (SPDC) [in is a light amplification process that takes place in a non- 
linear medium, where a photon from a pump field is converted into two photons, usually called signal and idler, with 
energy and momentum conservation. The signal and idler fields are therefore strongly correlated in energy, emission 
times, polarization and momentum. SPDC is a very convenient tool to produce entangled states of photons, which 
have been used to test the foundation of quantum physics and are at the heart of quantum information processing and 
communications (see Q for a review). In the basic setup, the pump field is cw, and the output state is the so-called 
two mode squeezed state (see e.g. Eq. (5.64)). When the pump intensity is low enough, the output state is well 
described by a large vacuum component plus a two-photon state, a photon pair. Recently, physicist have started to go 
beyond this basic configuration. On the one hand, when the classical pump field is no longer cw but pulsed. In this 
case, the down-conversion process can take place only when a pump pulse is "inside" the crystal, thus providing an 
information about the time at which the down-converted photons are emitted. Of course, as a counterpart, coherence 
is lost in the frequency domain, since a pulsedpump field is not monochromatic. Non-trivial effects of a pulsed 
pump have been predicted Isf and observed ^ for photon pairs. On the other hand, more efficient sources and 
large pump intensities allow to produce an output field where the four -and more- photons components are longer 
negligible. The four-photon component of the field is of interest in quantum optics Q and in quantum information, 
since four-qubit entanglement can be obtained |MI^ITo|. But this component can be a nuisance as well, for instance 
in quantum teleportation, because its presence decreases the fidelity of the two-qubit Bell state measurement [Tlj or 
in two-photon interference experiments where it limits the visibility |l2{ |. 

In this paper, we investigate both experimentally and theoretically the physics of the four-photon component 
produced in down-conversion, with a pulsed pump field. We start with a qualitative description of what is to be 
expected. The two meaningful quantities are the duration of the pump pulse At and the coherence time of the 
observed down-converted photons t'^'\ The characteristics of the four-photon state are captured by the relation 
between P2 and P4, the probabilities of creating 2, resp. 4 photons. 

Let us consider first the limit tP'^ ^ At. A large number of independent SPDC processes can take place inside 
a pump pulse .13J. In this limit, any 2n photon state can be satisfactorily described as n independent pairs. The 
probability of creating n pairs in a given pulse is described according to a Poisson distribution of mean value fi 

P„ = e '^fj,"'/nl, as shown in Section |II| In particular, for small /i we have P4 ~ and the four-photon state 
corresponds to two independent pairs, labelled |2 EPR). 

The other limit, t"^^ 3> At, can be achieved by the use of femtosecond pump pulses and narrow filtering of signal 
and idler photons. This condition is mandatory for all experiments where photons created in different SPDC events 
must interfere at a beam splitter, in order to preserve the temporal indistinguishability 14] . In this case, the emission 
of a "second pair" is stimulated by the presence of the first one [l^ leading to an entangled four-photon state |4 entg) 
that cannot be described as two independent pairs. Because of stimulated emission, we have P4 = ■ 

In the present paper, we study the transition between these two extreme situations. We define a parameter x G [Oj 1] 
that allows to interpolate between the Poisson distribution and the statistics arising from stimulated emission according 
to 
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In Section ini we present an intuitive description of the physics in the language of quantum states. In Section ITTll we 
present a simple experiment that allows to measure x- In Section ll VI we give the full quantum-optical formalism to 
describe the four-photon component of the field, and work out an approximate solution which agrees well with the 
experimental data. In particular, we find that x depends on t^.^ and At only through their ratio 

r = if/At. (2) 



II. THE FOUR-PHOTON STATE 



A coherent down-conversion process produces the two-mode squeezed state 
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where T = tanh^ and C = cosh<f and ^ is proportional to the amphtude of the pump field (see Eq. (5.64)). The 
state \n; n) describes the field with n photons in the signal mode and n photons in the idler mode; in particular, 
|2; 2) = |4 entg) is the four-photon state described in the introduction. 

In the limit where we study photons whose coherence is much larger than the time-bin (the width of the pump 
pulse), that is t^^ ^ At, there is a unique coherent down-conversion process taking place in the crystal for each pump 
pulse, and in this case the four-photon component of the field is |4 entg). We expect x = 1- 

If tP'' < At, a number N « At/tP'' = 1/r of independent SPDC processes can take place inside a pump pulse. To 
simplify this discussion, we consider N as an integer. The state of the field reads then 

l*)"""^ = ^ E ^" |2^) = ^ + ^|2) + T'l'i) + ...) , (4) 

n 

where the |2n) is the un-normalized superposition of all the states containing 2n photons. Specifically: 

The two-photon component |2) is the sum of the N states that describe "one pair in process j and no pairs in the 
other processes", that is |2) = |1; 1) (g) |0)... |0) + ... + |0) |0)... (g) |1; 1). Since all the components are mutually 
orthogonal, we have (2|2) — N . The four-photon component |4) is the sum of two kind of terms: (i) The N components 
1 2; 2) (g) |0) • • • (g) |0) -|- |0) ® |0) • • • g) |2; 2), in which the second pair is created by stimulated emission; each of these 
gives rise to the correlations of |4 entg). (ii) The N{N — l)/2 components where one pair is produced in process j and 
another pair is produced in a different process j'. Each of these gives rise to the correlations of |2 EPR). Therefore 
(4|4) = N{N + l)/2, and by normalizing this component we can say that the "four-photon state" is 



l^^^l^)) = \l |4 entg) + ^ ^— j |2 EPR) . (5) 

Referring back to we can compute the probabihties of having two or four photon: P2{N) = Pi{N) = 

_i_ rp4 N(N^+i) ^ limit of very large number of independent processes N, the usual argument leads to the 

Poisson distribution Now we have all the elements to compute x and relate it to the description of the four- 

photon component. For simplicity, we put C = 1. Then from (QJ we obtain x = 2P4{N) / P2{N)'^ — 1 that is 

1 



X = 

and we can re-write the four-photon state as 



^ = r for tP'^«At (6) 



I*4(X)> = V if^ |4 entg) + y |2 EPR) ■ (7) 

This provides an intuition on the link between x, the experimental parameter r and the entanglement in the four- 
photon state. 



III. THE EXPERIMENT 



A schematic of the experiment that measures x is shown in Fig. ^ A non-degenerate type I parametric down- 
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FIG. 1: Schematic of the experiment used to measure the stimulation parameter x- See text for details. 

converter is pumped by a pulsed laser. PDC modes a and h are then separated deterministically, via their different 
wavelengths. We ignore photons in mode 6; in mode a, we detect coincidence counts between the outcomes of a passive 
coupler. This coincidence measurement post-selects the events in which at least four photons have been produced 
in the down-conversion processes. The idea is to compare the events where four photons are created in the same 
pump pulse with the events where one pair is created in one pulse and another one in the next pulse. In the first 
case, we detect a coincidence in a time window centered at At^et = 0. The coincidence count rate i?o of this peak is 
proportional to \Pi- The factor i is the probability that the two photons exit different modes of the beam splitter. 
In the second case, we detect a coincidence with At^et — ^Tiaser (time between 2 laser pulses). We restrict ourself to 
the case when a photon created in pulse n is detected by detector Da,+ while photon created in pulse n-|- 1 is detected 
by detector Da-- The coincidence count rate Rside is thus proportional to j[P2{I)Y . Consequently, we have 
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Here is a description of the experimental setup. A mode-locked femto-second laser, generating Fourier transform 
limited pulses at 710 nm, is used to pump a lithium triborate (LBO) type I non linear crystal. The time between two 
subsequent laser pulses is l^Tiaser = 13ns. Collinear non-degenerate photon pairs at telecom wavelengths (1310 and 
1550 nm) are created by parametric down-conversion. The created photons are then coupled into an optical fiber and 
separated deterministically with a wavelength division multiplexer (WDM). We ignore the 1550 nm photons and we 
send the 1310 nm photons to a 50-50 fiber coupler. The two outputs of the coupler arc connected to photon detectors, 
labelled Da,+ and Da,- in Fig. ^ These detectors are Ge and InGaAs avalanche photodiodes (APD), respectively 
operating in Geiger mode. The Ge APD is used in passive quenching mode, while the InGaAs is used in the so called 
gated mode. In order to reduce the noise, the trigger for the InGaAs APD is given by a coincidence between the Ge 
APD and a 1 ns signal delivered simultaneously with each laser pulse. The signal from one detector serves as START 
for a Time to Amplitude Converter (TAG) , while the signal from the other one serves as STOP. We can thus measure 
the arrival times and see directly the effect of stimulated emission when the photons are created in a same pulse. An 
interference filter (IF) of different spectral width AA (5nm, lOnm, 40nm, FWHM) can be placed after the crystal, in 
order to change the coherence length of the down-converted photons. The coherence time (FWHM) is calculated from 
AA assuming a gaussian spectral transmission for the IF: t^^ = 0.44^^^ The gaussian transmission is a very 

good approximation for AA =5,10 nm, but it is less accurate for AA=40 nm. The calculated t^^ for AA=40 nm might 
therefore be underestimated. The pump pulse duration At can also be varied, and is measured after the crystal with 
an auto-correlator. Note that the pump pulses are no more completely Fourier transform limited after the crystal, 
due to chromatic dispersion in the optical path. With the different IF and the different pump pulse durations, it is 
thus possible to vary the ratio r = j td. Fig. [3 shows a typical TAG histogram for two different configurations 
(i.e. pump pulse duration and IF), leading to different values of r. For each value of r, one can directly measure the 
X by comparing the number of coincidence counts i?o in the central peak and Rside in the side peak. Fig. |3| shows 
the measured x, as a function of the ratio r. The theoretical prediction described later is 
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which by the way reproduces the predictions of Section ^] in the limiting cases (x = 1 for large r and x = ^ for 
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FIG. 2: TAG histogram for two different ratios r. The straight line corresponds to r=2.5 (measured x = 0.95) and the dotted 
line to r=0.2 (measured x = 0-3) The two curves are artificially shifted for clarity. 




FIG. 3: Measured x coefficient as a function of the ratio r. The open circles are experimental points. The full fine is a plot of 
formula @. 



small r) . The agreement of the data with this prediction is satisfactory — note in particular that there are no free 
parameters in the model. 
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IV. THEORETICAL DESCRIPTION 
A. State after down-conversion 

We want to describe the state of the e-m field produced as the output of a down-conversion process driven by a 
classical field. The scheme of the calculation is standard: the output state is 

1*) = e-«|o) c {i-tn~^n')\o), (10) 

where Ti- ^ J Hj{t)dt and Hj{t) is the Hamiltonian describing the down-conversion process in the interaction 
representation. This scheme has been applied in Refs. 0, Q to the case in which the pump field is not a continuous 
cw wave but a pulse of finite spatio-temporal extension. In these Refs, the calculation was limited to the two-photon 
term. Wang et al. studied the four-photon term, for degenerate type-I down-conversion. Here we consider colinear 
emission (propagation along z) of non-degenerate photons in a type-I crystal. 

The Fourier-transform of the pump field is written ■\/7pQ!(a;). Following the same steps as in Refs. 0,0 we find 
Ti = \fl {j^ + A) , where / is proportional to the intensity /p of the pump field, and where is the two-photon 
creation operator 

^ \ dcJi (ia;2 a(^i + CJ2) *&(wi, W2) a^(a;i)a^(a;2) ■ (11) 



Here appears the phase-matching function $(a;i, 0^2), peaked around k, fix and W2 ~ ^^2 with fJi 7^ ^1 |18| . In the 
following for conciseness we shall write 

V{i^\,^2) = a{i^\^ ^2)^{^\t^i) ■ (12) 

Now we must insert H. into H1U|I . Since MfS) = and ^y^^jO) oc |0), after removing the vacuum component that 
obviously does not contribute to the detection, one finds 

I*) ~ iVlA^\0) + ^A^^\0). (13) 

B. Probabilities and the observed x 

From 1^1') produced at the down-conversion, we can calculate the probabilities P2 = /(0|^^t|o) and P4 — 
^{OlA^A^ |0) of producing exactly two, respectively exactly four, photons. For this, one makes use of 

[a{Lj),a^{Lo')] = iSiu^-u') (14) 

and of a corollary of this commutation rule, that reads 

{0\a{uj")a{Lo"')a^uj)a^{uj')\0) 
= Siuj - u;")6iLu' - J") 

+(5(cj - - J') . (15) 

For Pi, the result is formally P2 ~ I J dujiduji [|7'(a;i, 1^2)!^ + "^(^^1, W2)7^*(a;2, wi)] . However, the second term of 
the sum is always zero because of the phase-matching condition (non-degenerate photons): in fact, the ranges of 
frequencies in which the first and the second argument of and consequently of V, lead to a non-zero contribution 
do not overlap. From now onwards, we simply consider that uji and UJ2 are different, that is [a(a;i), 0^(0)2)] = 0. Here 
is the point when our calculation differs from the one of Ref. . In conclusion, we have obtained 

P2 = I J da;ida;2\ViuJi,uj2)\^ = IJ2 . (16) 

To calculate P4, formula (|15|l is applied to the operators acting on both modes a and b, and one obtains 

P4 = [(J2)' + J4] (17) 
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where, writing du = dujidujbduj'iduj2. 



J duJ2duJ2 J duJiV*{uJi,Lu'2)'P{uJi,LU2] 



(18) 



One can verify that < J4 < ( J2)^ 0- By comparison with eq. ^ we find xo = Ji/i'hY] here, the suffix "0" means 
that this is the value of x the absence of any filtering after the SPDC. 

Now we must move and describe our experiment. Two approaches are possible. The "brute force" approach consists 
in effectively describing all the details of the experiment: write the pump pulse consisting of two well-separated pulses 
(this configuration gives similar experimental predictions as the one presented in section II, but is easier to compute), 
have mode a evolve through the coupler, and finally compute the coincidence rate at the detectors. This calculation 
is lengthy, although not devoid of interest for the theorist; we give it as an Appendix. The second approach is more 
clever: we know that the experiment measures x foi' ^ single pump pulse .9(1^-'), in the presence of an interference filter. 
We can then apply all that we have just done to find immediately 

X = J!l{j!f (19) 
where, writing G{x,y) = g{x + y)^{x,y) and F{lu) for the intensity profile of the filter, we have 

Ji = j duidU2F{LOl)\G{uJl,L02)\'' , (20) 

J! = j duiF{u^)F{uj[)G{ujuUJ2)G{uo[,uj'2)G*{uj[,0J2)G*{uji,uj'2) . (21) 
Obviously, in the Appendix we obtain the same result. 



C. Explicit estimate 



We have just found the general formulae that describe the quantity x measured in our experiment. In this paragraph, 
we solve explicitly (|2()|l and H21(l using some crude approximations; the final result will be formula Q for x, that has 
been shown to fit the experimental data. 

We make the following hypotheses: 
(i) The pump pulse is Fourier-transform limited: g{ijj) = (27rAp)"^^'' Q-it^-^p) /4Ap^ 
(n) The filter has a gaussian profile: F{uj) = {2tt1^f)^''^ ^-(^^n^)^ /2AI ^ 

(iii) The coherence time of the photons is uniquely determined by the width of the filter: i^'* — A^^. Therefore we 
replace $(1.^1, CJ2) by 1, that is, G{uji,uj2) by g{wi +^^-'2). 

The big advantage of this set of hypotheses however is that we are left with two meaningful quantities: Ap and 
Ap, whose inverses are the coherence times of the pump and of the photons in mode a. Plausible physical arguments 
will allow us to get rid of all these hypotheses at the end of the calculation. 

The calculation of jf is straightforward: via the change of variables (0^1,(^2) — > {^ii£, = ^1 +^2), the double 
integral factorizes into the product of two normalized gaussians, so Jf = 1. This implies x = J4 ■ The calculation of 
this integral is longer but not difficult: by the usual technique of square completion, one first arranges the terms in 
order to integrate out UJ2 and 0^2; the result of this procedure, writing lui — x and lu'i = y, is 



J[ = I dxdy exp 



2 V r2 a: 



with j4- = + 2S^- Using square completion again, one can integrate on y first and on x later. The result is 
jf" = ^ with = ]4- — A little more algebra leads to: 

From this simple result, we guess the general result Q by identifying f with r. This step is motivated by the following 
considerations. One the one hand, since we took a Fourier-transformed limited pulse for the calculation, the coherence 
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time of the pump 1/Ap is equal to the pulse duration, which we know to be the relevant quantity T^l. On the other 
hand, 1/AF is the coherence time of the down-converted photons, as long as the filter bandwidth is much smaller than 
the bandwidth of the unfiltered photons. When this condition is no longer satisfied, the relevant physical parameter 
is of course t^'' . 



V. CONCLUSION 



We studied the physics of four-photon states in pulsed parametric down-conversion. The parameters of the 
experiment determine to which extent the four-photon state exhibits four-photon entanglement, or can be rather 
described as two independent pairs. Any intermediate situation is quantified by a single parameter x that depends 
only on the ratio between the coherence time of the created photons and the duration of the pump pulse. A simple 
experiment to measure x has been realized. A theoretical model based on the standard formalism of quantum optics 
has been derived, that fits well the experimental data. Beside its fundamental aspect, this work is of practical 
interest in quantum optics, because it provides a simple mean to quantify the "coherence" of four photons states 
which is important for experiments such as quantum teleportation where independent photons must overlap at a 
beam splitter. 

Note : For related independent works on the statistics of photon numbers in down-conversion, see 

The authors would like to thank Claudio Barreiro and Jean-Daniel Gautier for technical support. Financial support 
by the Swiss NCCR Quantum Photonics, and by the European project RamboQ is acknowledged. 



VI. APPENDIX 



In this appendix, we want to re-derive formulae H20() and (|21|l with a full calculation of the experiment sketched in 
Fig. ID 

For the experiment we are going to consider, the classical pump field will be composed by two identical pulses 
separated by a time r. This configuration leads to the same experimental results as the one presented before and 
is easier to compute. If g{uj) is the Fourier transform (FT) of each pulse, the FT of the pump field is then simply 
a{uj) = g{uj){l + e^'^'^). To avoid multiplying notations, we keep V as in (|12|l for this explicit form of the pump field: 

V{uj,Cj) = (l + e*('^+'^)^)5(w + tD)$(w,tD) 

^ (l + e'("+-)-)G(w,w). (24) 



A. Evolution 



As discussed, the photons are separated according to their frequency. Those whose frequency is close to Q.i (resp. 
0,2) are coupled into the spatial mode a (resp. h). Photons in mode h do not undergo any evolution, while mode a 
evolves through a 50-50 coupler into modes c = (a, +) and d = (a, — ) according to 

a'(w) — * • (25) 

Inserting this evolution into l|13|l . the state at the detection reads 

m,,, ^ ^m + \Je\Q) (26) 



where, omitting the multiplicative constant I /2, we have written 

IC = [ dujdCjV{uj,Cj) [c\uj) + id\uj)\ b\Cj) . 



(27) 
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B. Detection (I): generalities 



We turn now to the detection. The experiment that we are describing involves the detection of two-photon coinci- 
dences. Let Dc and be the two detectors that monitor modes c and d. The probabihty of detecting a coincidence 
between the two events "photon detected in Dj at the time tj", j = c, d, is given by 

PD.,DAtc,t,) = \\Ei+Ht,)El;\t,)\^),J\^ (28) 

where 

Ei+HQ = j duU{v)c{v)e-^'''^ 

E],^\td) - j dvfd{v)d{v)e-''''^ (29) 

are the positive frequency component of the electric field at time tj in the mode detected by detector Dj j weighted 
by the amplitude of the filter fj{v) put in front of each detector. From now on, as in our experiment where the filter 
was actually put before the WDM, we consider fc = fd = f\ the intensity shape of the filter is F{uo) = |/(ti;)p. 

Now we insert ()26f) into H28|) . As expected, the term linear in K, gives no contribution: if there is just one photon 
in modes c or d, no coincidence count can be obtained. Similarly, no contribution comes from the terms of the 
form {uj)c^ {uj') and (i^(aj)d^(w') in K? ^ where both photons are found in the same mode after the coupler. In the 
calculation of the non-zero terms, we systematically omit multiplicative constants from now on. Using 1)14(1 for modes 
c and d, one finds: 

eI+\Qe[+\u)\^),^, = I d^e-(-*=+-''^) J{uo) f[u')V{u:,^)V{u\^')b\Co)b\u'm , (30) 



where we have written dui = dui du ' duj duj ' . The probability PDc,Dd{ic,td) is the square modulus of this expression. 
Using H15|l for mode 6, one finds 

PD^,DAtc,U)^ I d^e-'l("-")*=+("'-^')*-)l [K{u:,u:',uj,ib'-u,Oj') + K{uj,,o',u',Cj-u,Q')\ (31) 



where we have written dtj — duj duj' dCj dcu ' duj duj ' and 

K{lo,lo',lo,lo';u,Q') = f[uj)f{u;')f[Cj)f{Lo')V{uj,u)V{uj',u')V*[Cj,Cj)V*[Cj\Co'). (32) 

The detection rate (counts per pulse) is obtained by integrating PD^,Dditc, td) over the time-resolution of the detector 
2AT: 



R{T,,Td) = / dt^dtdPo^M^cM)- (33) 

The time-resolution must be longer than the coherence time of the photons (otherwise, the selected modes cannot 
be seen), and shorter than the spacing between the pulses, to allow the resolution in time-bins. As the result of this 
integration, R{Tc,Td) has the same form as PD^DditcTtd) given in (|31|l . via the replacements 



e 



-^{u.-^)t, ^ ATe-*("-")^'^sinc[(tj-c2')Ar] (34) 
g-»(a;'-i,')td ATe~'^'^'^^">'^-'smc[{uj' -u>')AT]. (35) 

Even if it seems redundant here, for subsequent ease, it is convenient to write down explicitly 

RiT,,Ta) - i?i(T„Td)+i?2(r„rd) (36) 

where we have defined, writing S{x) — Arsinc(a;AT'): 



Ri{Tc,Td) = y d^e-■'^('"-'^)^=+("'-^')^'*)lx(c^,c^^tJ,c2'';w,w')^(w-^:^)^(w'-^i'), (37) 
R2{Tc,Td) = /(ia;e-^[("~^)^-+("'-'^')^'')lX(c^,c^',ti',w;w,w')^('^-'^)^(^'-'^')- (38) 



9 



C. Detection (II): meaningful times 

We have just found an expression for R{Tc, Td). Now, recall that the first time-bin is defined by T = 0, the second 
time-bin is defined by T = r. Therefore, we expect R{Ti^,Tj) to be significantly different than only when and 
Td take the values or r. In particular, the counting rate in the central peak is Rc = R{0,0) + R{t,t), while 
Riat = ^(0, t) = R{t, 0) are the counting rates in each of the lateral peaks. We want to recover all these results 
out of our general formula. In addition, we shall have manageable expressions for i?^ and Riat , allowing a fit of the 
experimental data. 

Let us start with a qualitative argument, that is enough for our purposes. Recall the expression H24|l of V. The 
spectral width of g is larger than (equal to, for Fourier-transform limited pulses) which in turn is much larger 
than since we want two well-separated pulses. This means that g{Lu) is almost constant in a frequency range of 
width T^^. The phase-matching function is also constant over such ranges, because its typical width is the inverse of 
the coherence time of the down-converted photons. Now, suppose that Tc and Td are or r. If one inserts H24|l 
into the expression for R(Tc,Td) and develops the products, one finds that R{Tc,Td) is a sum of terms that are the 
product of g, <& and a phase factor of the form e*^'^, with ft some algebric sum of the w's. The arguments above prove 
that this phase fluctuates very rapidly in the frequency space, unless $7 ~ 0. Therefore, when we integrate over the 
co's, all the terms will average to but those whose phase factor is 1. Moreover, by direct check one can easily get 
convinced that if either of Tc or Td is equal to a time when no photon is expected, say ^ or 2t, then no phases can be 
erased: the coincidence rate becomes zero. In summary, the first step to simplify R{Tc,Td) consists in writing down 
explicitly all the terms, and keep only those terms whose phase factor is 1. From now onwards, we admit that Tc and 
Td are either or r. 

Having erased terms that fluctuate as r^^ in the frequency space, a further simplification is possible. The argument 
of the cardinal sine functions is {uj — (2>)Ar ^ AT/tP'\ But as we said, AT must be larger than t^'*, otherwise the 
photon cannot be seen by the detector. Therefore the cardinal sine will only be significant if w ~ (2), that is, we can 
replace sinc[(cL; — (2;) AT] with 6{uj — uj). 

All this simplification procedure is just a matter of patience. One finds that i?i(0, 0) = i?i(0, r) = Ri{t,0) = 
Ri{t,t) = {Jif given in while i?2(0,0) = R2{t,t) = J[ given in and i?2(0,T) = i?2(T,0) = 0. In 

conclusion, the detection rates in the central peak and in each of the lateral peaks are: 

Rc = i?(0, 0) + i?(T, r) - 2[( Jf )2 + J4] , (39) 
Riat = R{0,t) = R{t,0) = iJi)\ (40) 

and the ratio Rc/2Riat is equal to 1 + x as announced. 
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